Planar pyrochlore: a strong-coupling analysis 
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Recent investigations of the two-dimensional spin- 1/2 checkerboard lattice favor a val ence bond 
crystal with long range quadrumer order [J.-B. Fouet et at, preprint 3ond-mat/010807C |. Starting 
from the limit of isolated quadrumers, we perform a complementary analysis of the evolution of the 
spectrum as a function of the inter quadrumer coupling j using both, exact diagonalization (ED) 
and series expansion (SE) by continuous unitary transformation. We compute (i) the ground state 
energy, (ii) the elementary triplet excitations, (iii) singlet excitations on finite systems and find very 
good agreement between SE and ED. In the thermodynamic limit we find a ground state energy 
substantially lower than documented in the literature. The elementary triplet excitation is shown to 
be gapped and almost dispersionless, whereas the singlet sector contains strongly dispersive modes. 
Evidence is presented for the low energy singlet excitations in the spin gap in the vicinity of j = 1 
to result from a large downward renormalization of local high-energy states. 

PACS numbers: 75.10.Jm, 75.50.Ee, 75.40.-s 



Quantum-magnetism in low dimensions has received 
considerable attention recently due to the discovery of 
numerous materials with spin- 1/2 moments arranged in 
chain, ladder, and depleted planar structures. Many of 
these systems exhibit strongly gapped excitation spectra 
induced by various types of magnetic dimerization and 
discrete symmetry breakings. Yet, low energy collective 
spin dynamics may emerge from materials with strong ge- 
ometrical frustration leading to ground states with near 
macroscopicjiegeneracy and possibly even zero tempera- 
ture entropynEl. In this respect quantum spin systems on 
the kagome, and more recently on the pyrochlore lattice 
are of particular interest. While their classical counter- 
parts have been studied in considerable detail the role of 
quantum fluctuations in such systems is an open issue. 
For the spin- 1/2 kagome lattice gapless singlet excita- 
tions and a high density of singlet states in the singlet- 
triijJet gap have been established on finite systems by 
EEa. Similar analysis of the pyrochlore quantum-magnet 
is severely constrained by its three-dimensional stru^ 
Therefore, and as a first step, several investigations 
have focussed on the planar projection of the pyrochlore 
quantum-magnet, i.e. the spin- 1/2 checkerboard lattice 
of Fig. |l]a) for the case of = 1. 

In particular, ED at j = 1 has resulted in a sizeaiJe 
spin gap and a large-number of in-gap singlet statesBH. 
Moreover, recent EDO has given strong evidence in favor 
of a valence bond crystal (VBC) ground state with long 
range order in the 5* = quadrumers shown in Fig. [T^) . 
While this has been concluded from the quadrumer corre- 
lation at J = 1 it emphasizes the need for a perturbative 
investigation of the checkerboard magnet starting from 
the limit j ^ in Fig. ^). Low-order perturbation the- 
ory haaJaeen performed previously for the checkerboard 
magneto however with the unperturbed Hamiltoaian cho- 
sen at variance with the VBC found in the EDtl. There- 
fore, the purpose of this work is to explore the evolution 
of the spin spectrum of the checkerboard spin- 1/2 magnet 
as a function of the inter quadrumer coupling j, contrast- 
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FIG. 1: a) The checkerboard lattice. Spin- 1/2 moments are 
located on the open circles. Full (dashed) lines label the 
quadrumer bonds (bonds corresponding to the expansion pa- 
rameter j). b) Energy (E), spin (S), and local qi quantum- 
number of a single quadrumer. 



ing exact diagonalization (ED) against high-order series 
expansion (SE). 

Written in a form adapted to the VBC symmetry 
breaking and normalized to an overall unit of energy the 
Hamiltonian of the checkerboard magnet reads 
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where Pi...ji = S^i + . . . + Sji and S^i refers to spin- 
1/2 operators residing on the vertices i = 1 . . .4 of the 
quadrumer at site 1, c.f. Fig. |l|a). Hq is the sum over 
local quadrumer Hamiltonians the spectrum of which, c.f. 
Fig. [^b) , consists of four equidistant levels which can be 
labeled by spin S and the number of local energy quanta 
q\. Hq displays an equidistant ladder spectrum labeled by 
Q qi. The Q = sector is the unperturbed ground 
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state |0) oi Hq, which is the VBC of quadrumer-singlets. 
The Q = 1-sector contains local S ^ 1 single-particle 
excitations of the VBC with q\ ~ 1, where 1 runs over 
the lattice. At Q = 2 the spectrum of Ho has a total 
S' = 0, 1, or 2 and is of multiparticle nature. For S — 
at Q = 2 it consists of one-particle singlets with q\ = 2 
and two-particle singlets constructed from triplets with 
Qi — Qm = 1 and 1 ^ m. Consequently the perturbing 
terms oc j in (|l]) can be written as a sum of operators r„ 
which nonlocally create(destroy) n > (n < 0) quanta 
within the ladder spectrum oLiZo. 

It has been shown recentlju'ElB that problems of type 
allow for high-order SE using a continuous unitary 
transformation generated by the flow equation method 
of Wegneril The unitarily rotated effective Hamiltonian 
Hcs readsucl 
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where m = (toi . . . m„) is an 71 = |m|-tuple of integers, 
each in a range of G {0, ±1, . . . , ±A^}. In contrast to 
if of (0), H,s conserves the total number of quanta Q. 
This is evident from the constraint M{m) = Yl7=i ~ 
0. The amplitudes C(m) are rajLkmal numbers computed 
from the flow equation methodQtl. 

Previous applications of this method to spin models 
have been confined to dimer systemsEI where N < 2. Here 
we report on its first application using^a quadrumer ba- 
sis where N < 6. Explicit tabulationllJ of the r„ shows 
that for the checkerboard magnet = 4. Because the 
quadrumer basis consists of 16(4) states(spins) rather 
than only 4(2) states(spins) as for dimer systems, the 
maximum order achievable by a quadrumer expansion is 
lower in general than for a dimer expansion. Here we 
present results up to 0{7) for the ground state energy 
and the elementary triplet as well as up to 0{6) for the 
Q = 2 singlets. The C(m)-table is available on requestliil. 

Q-conservation leads to a ground state energy of 
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Evaluating this matrix element on clusters with periodic 
boundary conditions, sufficiently large not to allow for 
wrap around at graph-length n one can obtain SEs valid 
to 0{n) in the thermodynamic limit, i.e. for systems of 
infinite size. 

Q-number conservation guarantees the Q = 1-triplets 
to remain genuine one-particle states. A priori single- 
particle states from sectors with Q > 1 will not only dis- 
perse via iicffj but can decay into multi-particle states. A 
•posteriori however such decay may happen only at high 
order in j leaving the excitations almost true one-particle 
states. The dispersion of the single-particle excitations is 
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where t^^im = {fJ.Jm\Hcs\fJ,, 00) - SiynfloEf" are hop- 
ping matrix elements from site (0,0) to site {l,rn) for 
a quadrumer excitation fi inserted into the unperturbed 
ground state. For the thermodynamic limit t^^im has to 
be evaluated on clusters with open boundary conditions 
large enough to embed all linked paths of length n con- 
necting sites (0,0) to {l,rn) at 0{n) of the perturbation. 
(0|iioff|0) on the t^i^oo-chister. 
To assess the convergence of the SEs, as well as to 
identify the remaining excitation spectrum, we will con- 
trast the SEs with ED of finite systems. For this purpose 
eqns. (j^^ apply equally well, with however the cluster 
topologies set by the finite systems. We emphasize, that 
for some of the SEs generated by (||) on finite systems 
we were also able to compute standard SEs using a com- 
pletely different codala. In all cases we found both SEs 
to agree which serves as an independent check. 

Fig. 1^ summarizes our results for the ground state en- 
ergy per spin eg{j), the spin gap ET{.j), and the spectrum 
of the low lying singlets. The figure compares the SEs 
with our ED data on 2(3) x 2-4 systems, i.e. 16|(24) spins, 
as a function of j, as well as available ED dataS on 3 x 3 • 4 
sites, i.e. 36 spins, at j = 1. SEs are shown both for the 
identical system sizes and also for the thermodynamic 
limit to 0{j'^) for eg and Et as well as to O(j^) for the 
Q — 2-singlets. All SE data refers to the actual series, 
no Fade continuations have been applied. 

First we focus on the ground state energy in Fig. |^). 
In the thermodynamic limit we find 
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which is lower than the decoupled quadrumer-value of 
eg(0) = -l/2foraUj < 1 and yields e^(l) « -0.513677. 
This is also lower than the value of e^(l) w -0.49622 
which is above the bare quadrumer limit and has been 
reported from a difEerent third-order SE starting from 
decoupled tetrahedrm. While the SEs are not related to 
a variational principle the preceding is consistent with the 
decoupled quadrumers to incorporate the proper ground 
state correlations for j < 1 rather than the tetrahedral 
limit. This is confirmed on the finite systems in Fig. ^), 
where practically quantitative agreement is found for 
j < 0.7 between ED and SE for 2(3) x 2 - 4 sites. At j ~ 1 
very small deviations occur, which however decrease with 
system size. They are ^ 2% on the 16(24)jspin lattices 
and only 0.7% on the 3x3-4 site systemO. Regarding 
the SEs this agreement is remarkable also, since the high- 
order decrease of SE-coefficientalil tends to be slowed 
down on finite geometries. On the 2(3) x 2-4 lattice 
we observe a kink in the eg{j) as obtained from the ED 
at j « 1 which suggests a change of the nature of the 
ground state on these systems in the vicinity of j = 1 . 
This ispprobably a special property of the 16(24)-spin 
latticeda. 

Next we discuss the elementary triplet excitations. 
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36 spins. A clear cusp is observable in ET{j) close to j = 
1 on the 2(3) x 2-4 systems. This is closely related to the 
previously anticipated change in the nature of the ground 
state on the 16(24)-spin systems at j « 1 and is not 
captured by the SE. While no j-scan of Et is available 
on 3 X 3-4 sites it is remarkable that at j = 1 the SE agrees 
better with ED for 36 spins than for 16(24) spins and that 
the SE behaves qualitatively different, both on the 3x3-4 
lattice and in the thermodynamic limit as compared to 
the 2(3) X 2-4 lattices. I.e., on the larger systems the gap 
decreases monotonously, while it increases beyond j ^ 1 
on the smaller two. This suggests an absence - or strong 
reduction - of the discontinuities close to j = 1 on larger 
systems improving the convergence of the SE. 

Finally we discuss the low energy singlet excitations. 
At j -C 1 they result from the Q = 2 sector which involves 
solving a two-particle problem coupled to a one-particle 
problem. We have found the mixing between the one- 
and two-particle states in the S — 0, Q = 2 sector to 
be very weak. In fact, on the 2x2-4 system the mix- 
ing vanishes identically in the ED as well as the SE up 
to 0{j'^). For all other systems analyzed by SE includ- 
ing the thermodynamic limit, the mixing sets in only at 
0{j^). Therefore the Q = 2 one-particle singlet remains 
an almost true one-particle state even at finite j. In 
addition its dispersion E'^^^^ up to 0{j'^) is determined 
exactly by (^. In the thermodynamic limit we get 



FIG. 2: a) Ground state energy per spin, b) spin gap and 
c) singlet excitations with k = as a function of j. Lines 
correspond to the series, symbols to exact diagonalization. 
The numerical data at j = 1 for 36 spins is taken from Ref. ^ 
The grey shaded area in panel c) indicates that our numerical 
data is not complete in this region. 



The SE to 0{j'^) in the thermodynamic limit for the 
triplet dispersion reads 
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Remarkably, the hopping amplitude is exceedingly small 
and sets in only at 0(/) hinting at a fairly extended 
polarization cloud necessary in order to allow for triplet 
motionllj. From (Q) the spin gap occpjs at k = 0. 

The spin gap is shown in Fig. |^)Efl. Again, on finite 
systems of 2(3) x 2 - 4 sites and for j < 0.7 the agreement 
between ED and 0(j ^)-SE is very good. At j = 1 and for 
the SEs shown the absolute value of the relative difference 
to the numerically exact values for the spin gap improves 
from 11%, to 8.6%, to 2.7% on passing from 16, to 24, to 



In sharp contrast to (|^) this low energy singlet is strongly 
dispersive. Its gap resides at k = 0. 

Fig. &;) shows the evolution of all eigenvalues of H 
from (pp^ in the singlet sector on the 3x2-4 system at 
k = as well as the spectrum of iJcff at 0(/) in the 
5 = 0, Q — 2 sector. For j < 0.7 the agreement between 
ED and SE is satisfying. Similar agreement is found at 
the remaining k- vectors. The lowest S = excitation can 
be identified with a singly degenerate Q = 2 two-triplet 
excitation. The binding energy for this two-particle state 
is approximately proportional to j with a coefficient of 
order unity. Starting at j ^ 0.7 a large number of singlets 
from higher Q-sector enter the low energy spectrum. At 
j — 1 several of these singlets reside in the spin gap. This 
is consistent with ED at j = 1 in Refs. The singlet 
gap occurs at k = 0. Preliminary analysis indicates that 
the lowest singlets have a substantial overlap with the 
(5 = 4 sector. The numerical spectrum in Fig. ^) is 
very suggestive of a softening in the singlet sector in the 
vicinity of j = 1. If this is a feature particular only to the 
finite systems is an issue in clear need of further study. 
Close to j = 0.7 a near level crossing occurs, with the 
lowest state of the SE for j > 0.7 corresponding to a 
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FIG. 3: System size dependence of the dispersion from SE 
for Q = 1, 5 = 1 and Q = 2, S = 0. Lines: ErCkJ) (al- 
most nondispersive) and i?s,ipt(k, j) (dispersive). Symbols: 
all eigenvalues of the Q — 2, S = sector. 

continuation of the Q — 2 one-particle singlet. 

Finally we speculate on the finite size dependence to 
be expected for systems presently inaccessible to ED fo- 
cusing on the Q = 1 triplet- and the Q ~ 2, S = 
singlet-sector by using SE instead. In Fig. | we show 
the spectrum of i?eff for these two sectors at O(j^) for 
j — 0.8 and 1.0 along a path in the Brillouin zone (BZ) 
on four systems ranging from = 3 x 3 • 4 sites up to 
the thermodynamic limit. As we have pointed out SEs 
will converge more rapidly on larger systems. Therefore, 
at j = 0.8 the SE resuhs of Fig. | will be rather well 
converged while at j = 1 they will be qualitatively cor- 
rect. Fig. ^ includes the one-particle triplet and singlet 
dispersions, i?T(k,j) and i?s_ipt(k, j) with k varying con- 
tinuously along the path in the BZ. For the latter curve 
Es,iptO^, j) has been evaluated at 0{j^). Moreover in the 
thermodynamic limit and for Q — 2 only Es^ipt()<-, j) is 
depicted. First Fig. |^ illustrates the difference in strength 
of the dispersion between the elementary triplet and al- 
most one-particle singlet from the Q = 2-sector. Sec- 



ond the figure shows that both, the elementary triplets 
and the corresponding, almost two-particle singlets with 
Q = 2 display only a small downward shift with increas- 
ing system size. However, the spectrum of the almost 
one-particle singlet with Q = 2 experiences a strong re- 
duction of its band-width including an upward shift of 
its gap. In particular, at j = 1 the SE suggests that the 
Q — 2-singlets are above the spin gap for sufficiently large 
systems and that only singli^s from sectors with Q > 2 
might reside in the spin gapEj. This underlines the need, 
both, for larger systems in ED and larger Q-values in SE 
studies of the singlet sector. 

In conclusion, by comparing ED and SE we find the de- 
coupled quadrumer-limit to be an excellent starting point 
for a perturbative treatment of the spin- 1/2 checkerboard 
magnet even at j w 1. The remarkably fast convergence 
of the series for the ground state and the elementary 
triplet excitation up to at least j = 1 indicates that the 
points j = and 1 belong to the same phase. This is 
consistent with the strong VBC correlations observedO at 
j = 1. We have further shown the elementary triplets of 
the checkerboard magnet to be almost nondispersive. Up 
to intermediate j the elementary singlet excitations are 
almost nondispersive two-particle states from the Q = 2 
sector residing above the spin gap. On finite systems and 
in the vicinity of j = 1 , we find a strong downward shift of 
singlets from higher Q sectors, entering the spin gap. In 
future work it would be interesting to study these higher 
Q-sectors by SE and to obtain j-dependent ED spectra 
on larger systems, extending in particular the 36-site re- 
sults of Ref. H to j 1. 
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We have performed a similar analysis for all distinct k- 
sectors on the 2(3) x 2 ■ 4 systems. We note that on the 
2x2-4 lattice i5j^'^"'(k, j) is completely k- independent, 
both, in ED and SE to 0{f). 

Since the jVBC picture implies a two-fold ground state 
degeneracyO at j = 1, we expect at least one low- lying sin- 
glet excitation to come from these higher Q-sectors. Due to 
Q-conservation it is possible that the mixing of this state 
with the ground state as well as the convergence problems 
usually associated with this can be avoided in the present 
case. 



